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A.  Introduction 

Polymers  composites  are  used  as  a  protection  against  ballistic  impact  for  the  buildings  and 
military  vehicles.  The  main  function  of  an  armor  is  to  make  the  impact  of  the  projectile  to 
cause  no  damage  to  the  functionality  of  the  structure.  Initially  when  the  armors  were  made  this 
function  was  performed  by  putting  a  large  mass  on  all  the  vulnerable  regions.  It  works  on  the 
principle  of  conservation  of  momentum.  The  buildings  were  protected  by  using  cheap  material 
with  high  thickness.  The  drawback  of  this  approach  is  that  in  order  to  meet  higher  intensity 
blasts  more  material  must  be  used  and  this  resulted  in  further  increase  in  the  weight.  This 
results  in  more  difficulty  in  the  transportation  of  these  material  and  reduction  in  the  speed 
of  the  vehicles.  When  designing  armor  weight  started  to  play  a  important  role  addded  to 
functionality.  Hence  the  need  for  the  use  of  the  lightweight  structures  to  perform  the  function 
arose.  Polymer  composites  satisfied  both  the  functionality  and  weight  requirements.  Nowadays 
these  polymers  are  also  used  as  coating  for  the  building  instead  of  the  conventional  thicker  wall 
for  the  ballistic  protection.  There  commercial  polymers  such  as  PAXCON™,  Dragonshield- 
HT™,  Rhino  Armor™  which  use  polyurea(PU),  polycarbonate(PC),  polymethylmethacrylate 
(PMMA)  etc.  for  the  ballistic  protection.  For  example,  it  can  be  seen  from  the  figures  1 
through  3  that  spraying  a  coating  of  RhinoArmor  protected  the  wall  from  the  blast.  From 
figure  4  it  can  be  seen  that  Dragonshield-HT  has  been  used  as  protective  coating  for  military 
vehicles. 

The  hypothesis  is  that  it  is  possible  to  use  a  two  network  theory  to  capture  the  essential 
features  of  a  impact  response  of  multilayered  polymers  and  solve  the  resulting  boundary  value 
problem  using  finite  volume  scheme. 

Traditionally  in  the  explicit  constitutive  theory  for  the  viscoelastic  solids,  the  current  stress 
is  a  functional  of  the  strain  history.  As  discussed  by  Rajagopal  and  Srinivasa  [1]  the  use  of  this 
theory  results  in  more  difficulties  in  estimating  the  model  parameters  since  it  is  impossible  to 


2 


Fig.  1.  Wall  on  the  left  is  coated  (Rhino Armor) 
and  wall  on  right  not  coated 


Fig.  2.  Blast  made  in  front  of  the  wall 


control  the  different  components  of  strain  and  also  most  of  the  experiments  are  either  partially 
or  fully  stress  controlled.  While,  in  the  implicit  constitutive  theory  there  exists  a  functional 
relationship  between  the  stress  and  strain.  This  representation  has  the  following  advantages 

•  The  response  functions  of  the  model  can  be  easily  compared  against  experimental  data 
and  model  parameters  obtained. 


•  It  results  in  simple  ODE  which  is  easily  amenable  to  numerical  scheme. 

Hence,  in  this  study  we  use  a  implicit  constitutive  theory  for  the  mechanical  response  of  a 
viscoelastic  solids. 


B.  Advantages  of  using  multilayered  plates 

Polymers  such  as  PMMA,  PC,  PU  etc.  are  layered  together  to  improve  the  performance  of 
structures  which  are  subjected  to  impact  loading[2,  3,  4,  5].  Many  experimental  and  numerical 
studies  have  been  carried  out  to  study  the  response  of  polymers.  For  example  experiments  have 
shown  that  the  performance  of  a  bilayer  plate  under  impact  was  improved  by  the  presence  of 
a  polymeric  material[5].  The  performance  depends  not  only  upon  the  type  of  the  polymeric 
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Fig.  3.  The  coated  wall  does  not  suffer  any  se¬ 
rious  damage  while  the  uncoated  wall  is 
broken  completely  exposing  the  interior. 


Fig.  4.  Polyurea  coating  to  protect  military  ve¬ 
hicle 


material  but  also  upon  the  thicknesses  and  layering  sequence.  For  example,  bilayer  experimen¬ 
tal  results  in  [2]  suggests  that  the  presence  of  a  polymeric  material  on  the  impact  side  reduces 
the  performance  when  compared  to  its  presence  on  the  other  side.  The  presence  of  a  com¬ 
pliant  material  between  two  stiff  layers  increases  the  performance  by  about  100%  compared 
to  the  individual  layers  whereas  a  bilayer  composite  increases  the  performance  by  about  25% 
[5].  Sudden  impact  of  these  material  results  in  the  generation  of  waves  and  these  waves  gets 
attenuated,  reflected  and  dispersed.  The  speed  at  which  the  wave  moves  depends  upon  the 
type  of  the  material.  Wave  propagation  through  a  homogenous  material  and  layered  composite 
are  explained  by  considering  a  single  wave(figures  5(a),  5(b)  and  5(c)).  In  a  homogenous 
material,  the  wave  propagates  in  the  same  direction  without  any  reflections  but  gets  attenuated 
(figure  5(a))  and  dispersed.  In  a  layered  composite,  as  the  wave  hits  an  interface  2,  it  partly 
gets  reflected  and  partly  gets  transmitted.  For  example  if  we  assume  that  the  PMMA  has  a 
lower  elastic  modulus  than  PC,  then  the  velocity  of  the  transmitted  wave  is  greater  in  PC  (high 
elastic  modulus)  when  compared  to  that  of  PMMA(low  elastic  modulus)  which  is  represented 
schematically  in  figure  5(b)  and  5(c).  In  order  to  study  phenomena  such  as  this,  as  well  as  to 
develop  a  tool  that  will  enable  us  to  quickly  evaluate  the  impact  resistance  of  different  kinds 
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of  polymers  layered  in  a  variety  of  ways,  it  is  necessary  to  develop  a  physically  realistic  and 
computationally  efficient  model. 
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the  wave  trav¬ 
els  without  any 
reflection 
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reflected  waves  at  the 
interface  2  at  different 
times.  The  former 
travels  at  higher  speed 
in  PC  compared  to  the 
latter  in  PMMA 
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to  the  latter  in  PC 


Fig.  5.  Schematic  representation  of  the  propagation  of  a  single  wavefront  in  different  material 
compositions.  The  location  of  the  wavefront  at  the  interfaces  at  three  different  times  are 
shown  by  the  arrows. 


C.  Numerical  scheme 


The  Finite  volume  method(FVM)  was  used  in  this  study  to  numerically  solve  the  govern¬ 
ing  equations.  In  this  procedure,  the  whole  domain  is  divided  into  a  finite  number  of  non¬ 
overlapping  control  volumes.  For  each  control  volume,  the  mass,  momentum  and  energy 
balance  equations  are  exactly  met.  The  constitutive  equations  are  approximated  by  suitable 
interpolation  schemes  across  neighbors.  With  this  approach,  it  is  possible  to  meet  local  conser¬ 
vation  of  mass,  energy  and  momentum  and  only  the  response  functions  are  approximated.  In 
the  finite  element  method(FEM)  the  local  conseri’ation  is  met  only  when  the  weighting  function 
is  independent  of  time  or  when  using  a  discontinous  Galerkin  method[6].  The  finite  volume 
method(FVM)  has  been  very  successful  in  computational  fluid  dynamics(CFD)  but  it  has  not 
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been  extensively  used  like  FEM  in  the  area  of  solid  mechanics.  Studies  have  been  carried 
out  to  predict  the  linear  elastic  behaviour  of  a  incompressible  elastic  solid  by  using  the  finite 
volume  method  [7,  8].  This  procedure  was  used  to  discretise  the  momentum  balance  and  kine¬ 
matic  condition  of  incompressibility  and  hence  the  three  dimensional  small  strain  problem  was 
solved  using  an  algorithm  similar  to  SIMPLE™ [9].  The  FVM  was  also  used  to  calculate  the 
stresses  and  displacement  by  Wheel  [10]  for  the  axisymmetric  bodies.  The  stress  equilibrium 
equations  in  terms  displacement  gradient  was  obtained  by  using  the  linear  elastic  constitutive 
equations  and  linear  interpolations  was  used  for  the  gradients.  This  procedure  was  also  ex¬ 
tended  to  solve  small  strain  problems  in  linear  viscoelastic  solids  to  be  used  in  applications 
involving  heating  or  cooling  of  viscoelastic  materials.  The  integral  form  of  constitutive  rela¬ 
tion  between  stress  and  strain  histories  along  with  thermal  effect  was  assumed  and  FVM  was 
used  to  solve  the  momentum  and  energy  balance  equations[ll].  Slone  et  al.  [12]  studied  the 
dynamical  structural  response  of  linear  elastic  solids  using  three  dimensional  FVM  and  an  im¬ 
plicit  scheme  for  the  time  step.  The  discretisation  method  was  similar  to  the  previous  studies. 
Studies  discussed  above  used  FVM  to  solve  either  linear  elastic  or  linear  viscoelastic  solid 
problem.  As  discussed  by  the  authors  in  the  previous  studies[7,  10,  11]  FVM  has  the  following 
advantages 

•  simple  and  efficient 

•  method  is  conservative  on  local  and  global  scale 

•  simple  to  prescribe  the  boundary  conditions 


Hence  in  this  study,  we  used  FVM  to  solve  the  impact  response  of  layered  viscoelastic  solids. 
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D.  Modelling  approach 

The  constitutive  relations  for  viscoelastic  materials  are  usually  represented  in  two  different 
forms  -  integral  and  differential/rate  type[13].  In  a  typical  integral  model,  the  current  value  of 
the  stress  is  given  in  terms  of  the  integral  of  the  history  of  the  strain  rate,  with  a  kemel(  usually 
called  the  relaxation  modulus)  that  encodes  the  fading  memory  of  the  material.  Simulating 
such  a  model  is  a  challenge  ([14])  because  of  the  necessity  of  carrying  the  information  about  the 
history  of  the  material  response  over  time.  In  order  to  circumvent  this,  the  relaxation  moduli  are 
usually  represented  as  a  sum  of  series  of  exponentials  which  is  referred  to  as  a  Prony  series([15, 
16]).  The  Prony  series  have  many  parameters  that  are  difficult  to  determine  experimentally.  In 
the  case  of  non-linear  integral  models,  the  challenges  are  further  exacerbated. 

An  alternative  approach  is  the  use  of  differential/rate  type.  The  distinction  between  a 
differential  and  rate  type  model  is  based  on  the  representation  for  the  stress.  The  stress  in  the 
differential  type  model  is  defined  explicitly  in  terms  of  appropriate  kinematical  variables  and 
their  derivatives.  For  the  stress  in  rate  type  equations  there  is  a  relationship  between  the  stress 
and  its  time  derivatives  and  the  kinematical  variable  and  its  time  derivatives.  Thus,  fluids  of 
grade  n  (see  [17])  are  referred  to  as  differential  type  models  while  the  Oldroyd-B  fluid  (see 
[18])  is  referred  to  as  a  rate  type  model. 

E.  Wave  Propagation 

The  equation  of  propagation  of  waves  through  a  linearized  elastic  solid  in  one  dimension  is 
represented  as 

d2u  _  d2u 
^  dt 2  dx 2 


(0.1) 
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where  u,  p,  E,  x  and  t  represent  the  displacement,  density,  elastic  modulus,  space  and  time 
domain  respectively.  D’Alembert  solution  to  the  above  equation  is  of  the  form 


u(x,  t )  =  f(x  —  ct)  +  g(x  +  ct) 


(0.2) 


where  c  is  velocity  of  propagation  of  the  wavcf  y/ E/p).  A  similar  approach  was  also  used  for 
viscoelastic  solids  in  studies[3,  19,  20,  21,  22,  23,  24]  which  will  be  discussed. 

Most  of  the  models  considered  a  linear  viscoelastic  model  for  describing  wave  propagation 
in  the  polymers  beacause  obtaining  an  analytical  solution  would  be  comparatively  easier.  Lee 
and  kanter  [19]  arrived  at  the  solution  of  the  wave  propagation  in  a  Maxwell  fluid  by  comparing 
the  governing  equations  of  motion  with  the  telegraphic  equation  which  is  given  by 


1  d2u  _  d2u  (R  G\  du  GR 
LCdx2  ~W+\L  +  Cjdi  +  ~LCU 


(0.3) 


where  u,  R,  L,  G  and  C  are  volatage,  resistance,  inductance,  conductance  and  capacitance 
respectively.  The  equations  were  solved  using  the  Laplace  transform  technique  resulting  in  a 
solution  involving  Bessel  functions.  Theocaris  and  Papadopoulou  [20]  considered  the  strain 
to  be  harmonic  and  solved  the  one  dimensional  wave  equation  for  the  displacement  along  the 
axis  of  the  rod  where  the  solutions  involves  attenuation  and  wave  velocity  parameters.  The 
solution  of  the  wave  equation  in  this  case  involve  two  constants  representing  the  direction  of 
wave  propagation  i.e.  one  for  the  positive  direction  and  the  other  for  the  negative  direction.  The 
two  constants  depend  upon  the  frequency  of  the  propagating  wave.  Fourier  series  was  used  to 
numerically  calculate  these  constants.  The  impact  response  of  polymeric  materials  was  studied 
by  assuming  a  linear  viscoelastic  relation  for  the  polymeric  material  [21].  Musa  [22]  predicted 
the  impact  response  by  solving  the  governing  equations  of  a  linear  viscoelastic  material  using 
Laplace  transform.  The  scattering  and  dispersion  of  the  stress  waves  was  observed  experi¬ 
mentally  by  Zhuang  [23].  Experimental  and  numerical  studies  was  carried  out  by  Boteler  et 
al.  [3]  to  study  the  wave  propagation  in  a  polymer  matrix  composite  subjected  to  impact.  He 
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used  the  integral  model  and  the  Prony  series  to  express  the  kernel  was  obtained  from  Zocher  et 
al.  [24]  and  used  to  predict  the  experimentally  observed  phenomenon  of  wave  dispersion  and 
attenuation.  This  model  was  able  to  capture  the  attenuation  but  not  the  dispersion. 

Many  models  wherein  a  differential  equation  relates  the  strain  and  the  stress  can  be  ob¬ 
tained  within  the  context  of  thermodynamics  by  assuming  energy  storage  and  energy  dissipa¬ 
tion  mechanisms.  From  the  point  of  view  of  one-dimensional  modeling  these  energy  storing 
and  energy  dissipation  mechanisms  can  be  viewed  as  elastic  springs  and  viscous  dashpots, 
respectively.  Within  the  context  of  one  dimensional  nonlinear  response,  the  combinations  of 
linear  springs  and  dashpots  represent  the  energy  storage  and  energy  dissipative  mechanisms, 
respectively  [25].  In  the  development  of  three  dimensional  models,  one  dispenses  with  the  idea 
of  spring  and  dashpots  but  directly  uses  expression  for  the  same  from  which  the  response  rela¬ 
tions  can  be  derived  (see  [26,  27]).  Such  models  are  typically  purely  mechanical  in  nature  and 
do  not  take  into  account  thermomechanical  phenomena.  Rajagopal  and  Srinivasa  [1,  28]  have 
developed  a  three  dimensional  framework  based  on  the  maximisation  of  the  rate  of  entropy 
production  for  thermoviscoelastic  solids.  The  mechanical  response  of  polymers  is  modelled 
by  assuming  it  as  a  viscoelastic  solid  since  it  exhibits  an  instantaneous  elastic  response  for  a 
sudden  application  of  load  and  creep  behaviour  when  the  load  is  held  constant.  The  simplest 
model  to  capture  the  above  two  basic  response  is  two  network  theory  composed  of  an  elastic 
(backbone)  network  coupled  with  a  viscoelastic  transient  network  [1,  28].  Thus  the  total  stress 
a  is  the  sum  of  the  partial  stresses  <j\  and  a2  in  each  of  the  networks,  i.e., 


a  —  cti  +  cr2 


(0.4) 


This  theory  will  be  used  in  this  study  to  capture  the  mechanical  response  of  the  polymer  which 
are  subjected  to  impact. 
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F.  Two  dimensional  viscoelastic  model 


The  one  dimensional  model  is  discussed  in  [29].  Let  us  consider  a  plate  that  is  infinite  in  z 
direction  and  has  finite  dimensions  of  length  T  in  x  direction  and  thickness  'll  in  y  direction. 
From  the  use  of  of  implicit  contistitutive  framework  we  assume  that  the  response  of  material 
is  nonlinear  eventhough  the  strains  are  small.  The  cartesian  coordinate  system  was  used  in  this 
study.  The  stress  tensor  and  velocity  vector  at  any  time ’t’  are  represented  by  T(x,y,t),  v(x,y,t) 
repectively. 

The  staggered  grid  in  two  dimension  is  shown  in  figure  6  along  with  a  sample  control 
volume  where  ’T’  is  the  stress,  'v'  and  V  x  and  y  component  of  velocity.  From  conservation 
of  linear  momentum  we  know  that 

r\ 

P ^  =  divT  (0.5) 

The  total  stress  on  T  on  the  network  is  the  sum  of  the  partial  stresses  T1 1  j  and  T(2)  in 
each  networks,  i.e., 

T  =  T(1)  +  T(2)  (0.6) 


The  mechanical  response  of  network  1  is  given  by 

<otO)  _  dwk 
dt  ijkl  d x* 

The  mechanical  response  of  network  2  is  given  by 

^r  =  M^^_7?(T(2)*T(2))T(2) 


(0.7) 


(0.8) 


where  C ijki  and  are  fourth  order  elasticity  tensors. 
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Fig.  6.  Staggered  grid  in  two  dimension 
G.  Nondimensionalization 

Before  using  the  finite  volume  technique  to  solve  the  momentum  and  constitutive  equations 
it  is  necessary  to  choose  appropriate  scales  to  nondimensionalize  them.  The  main  reason  for 
the  nondimensionalization  is  to  reduce  the  number  of  independent  parameters  so  that  we  could 
assess  the  factors  influencing  the  solution.  The  non  dimensionalized  equations  used  here  are 
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H.  Solution  Technique 


The  commonly  used  numerical  techniques  to  solve  the  momentum  equations  are  finite  element 
and  finite  volume  methods.  Both  involve  the  integration  of  the  differential  equations  to  obtain 
the  weak  form.  The  main  advantage  of  FVM  over  FEM  is  conservation  of  quantities  at  each 
and  every  control  volume.  Integrating  the  x  momentum  equation  over  velocity  control  volume 
(see  figure  7) 


>■ 


r-e  w  rN  Qji 

xxdW  +  / 


dx 


dy 


(0.9) 


Now  y  momentum  equation  is  integrated  over  velocity  control  volume 


>  - 


a rp  pn  oth 

xydw  +  /  ?i™dw 


'W 


dx 


dy 


(0.10) 


Similarly  the  constitutive  equations  are  also  integrated  and  the  final  discretised  equations 
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Fig.  7.  Control  Volume  for  the  2d 
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I.  Material  parameters  for  the  model 

The  material  parameters  were  calculated  using  the  experimental  results  obtained  from  the  lit¬ 
erature.  PC  and  PU  were  considered  for  this  study.  PC  creep  compliance  curves  were  obtained 
from  [30].  PU  creep  compliance  curve  were  obtained  from  [31].  Creep  test  data  experiment 
for  PC  taken  from  literature  was  carried  out  at  a  constant  pressure  46.3  MPa  for  a  time  period 
of  3600  s  at  room  temperature.  For  PU,  the  experiments  were  performed  at  a  constant  pressure 
of  2.5  MPa  for  a  time  period  of  560  s  at  a  tern-  perature  of  30  C.  The  material  parameters  were 
calculated  by  minimizing  the  error  between  the  experimental  data  and  the  model  predicted  data 
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using  fminsearch  function  in  MATLAB.  The  parameters  obtained  are  shown  in  Table  1.  The 
fit  of  experiment  vs  model  for  PC  and  PU  are  shown  in  Figure  8  and  9  respectively. 

Table  I.  Material  Parameters  based  on  creep  compliance  data 


PC 

PU 

Ex 

1.7043  *  109Pa 

7.0658  *  107Pa 

e2 

3.8630  *  108Pa 

9.7342  *  107Pa 

V 

6.4173  *  10"  14 

6.4173*  10"  14 

5.6 
5.4 
5.2 

5 
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4.6 
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450  540 


Fig.  8.  Experiments  on  PC([30])  compared  with 
the  implicit  model  presented  here.  No¬ 
tice  the  excellent  match  for  the  nonlinear 
creep  compliance  based  on  the  constitu¬ 
tive  parameters  given  in  non-dimension- 
alized  momentum  equations. 


Fig.  9.  Experiments  on  PU([32])  compared  with 
the  implicit  model  presented  here.  No¬ 
tice  the  excellent  match  for  the  nonlinear 
creep  compliance  based  on  the  constitu¬ 
tive  parameters  given  in  non-dimension- 
alized  momentum  equations. 


J.  Results  and  Discussion 

In  the  study  conducted,  six  different  material  combinations:  (1)  pure  PC,  (2)  pure  PU,  (3) 
bilayer  (PU/PC),  (4)  bilayer  (PC/  PU),  (5)  trilayer  (PC/PU/PC)  and  (6)  trilayer  (PU/PC/  PU) 
were  considered.  One  face  of  the  plate  is  fixed  while  the  other  face  has  a  small  regions  which 
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is  given  a  ramp  stress  in  ’x’  direction  with  slope  ’a’  for  a  time  ’£0’  and  is  maintained  constant 
till  All  the  remaning  boundaries  are  stress  free. 

K.  Wave  Propagation 

Sudden  impact  of  a  material  results  in  the  generation  of  waves.  The  speed  at  which  these  waves 
move  depends  upon  the  type  of  the  material. 


L.  Stress  Propagation 


The  propagation  of  stress  at  different  instances  of  time  is  shown  in  Figures  10  through  15 
through  the  thickness  of  the  plate  when  impact  is  made  on  PC  and  PU.  Initially  the  stresses  are 
0  throughout  the  thickness  for  all  cases.  The  propagation  is  faster  in  PC  layer  when  compared 
to  that  of  PU  since  it  is  more  compliant  than  PU.  In  case  of  the  trilayer  multiple  reflections  can 
be  observed  from  the  interfaces. 


x  10-5 


Fig.  10.  Stress  propagation  at  20  units  of  time 
when  PC  is  impacted.  Transmission  of 
stress  through  homogenous  layer  PC 


x  direction  of  intial  stress  propagation 


Fig.  11.  Stress  propagation  at  30  units  of  time 
when  PC  is  impacted.  Reflection  and 
transmission  of  stress  through  homoge¬ 
nous  layer  PC 
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x  direction  of  intial  stress  propagation 

Fig.  12.  Stress  propagation  at  20  units  of  time 
when  PC  is  impacted.  Transmis¬ 
sion,  reflection  and  refraction  of  stress 
through  heterogenous  layer  PC/PU/PC 

M.  Stress  and  Kinetic  energy  on  wall 


x  direction  of  intial  stress  propagation 

Fig.  13.  Stress  propagation  at  30  units  of  time 
when  PC  is  impacted.  Multiple  reflec¬ 
tions  refraction  of  stress  from  wall  and 
interface  through  heterogenous  layer 
PC/PU/PC 


The  important  application  of  interest  with  regard  to  the  composite  polymer  that  is  used  here  is 
to  reduce  the  impact  on  the  wall  besides  which  it  is  in  contact.  The  peaks  in  Figure  16  and  17 
represents  the  time  when  the  impact  wave  reaches  the  wall.  The  damping  effect  on  the  wave 
can  be  seen  from  the  gradual  decrease  in  the  peak  response  as  the  wave  travels  back  and  forth. 
The  most  important  aspects  that  needs  to  be  taken  into  account  in  the  designing  of  polymer 
composites  as  shock  absorbers  is  the  magnitude  of  the  stress  and  kinetic  energy  on  the  the 
wall.  The  above  indicates  that  the  magnitude  of  the  stress  and  kinetic  energys  is  lowest  in  case 
of  the  trilayer  (PC/PU/PC)  with  regard  to  the  six  different  layers  that  were  considered. 
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x  10' 


x  direction  of  intial  stress  propagation 


Fig.  14.  Stress  propagation  at  20  units  of  time 
when  PU  is  impacted.  Transmis¬ 
sion,  reflection  and  refraction  of  stress 
through  heterogenous  layer  PU/PC/PU 

N.  Conclusion 


Fig.  15.  Stress  propagation  at  30  units  of  time 
when  PU  is  impacted.  Transmis¬ 
sion,  reflection  and  refraction  of  stress 
through  heterogenous  layer  PU/PC/PU 


The  performance  of  the  different  materials  such  as  PC  and  PU  under  impact  was  studied  us¬ 
ing  a  nonlinear  rate  type  model.  The  governing  equations  were  solved  numerically  using  a 
FVM  approach.  The  parameters  for  characterizing  the  model  were  obtained  from  the  creep 
compliance  data  found  in  the  literature.  The  current  one  dimensional  model  discussed  in  the 
paper  only  allows  for  transmission  and  reflection  of  the  waves.  The  one  dimensional  model 
discussed  in  the  paper  was  able  to  capture  the  reflection  phenomenon  due  to  the  interaction 
of  the  waves  at  the  interface  but  fails  to  capture  the  dispersion  phenomenon  which  the  two 
dimensional  model  was  able  to  capture.  The  comparison  of  the  stress  in  the  wall  when  using 
different  combinations  of  PC  and  PU  was  carried  out.  The  results  show  that  the  kinetic  energy 
was  lower  on  the  wall  when  using  a  composite  layer  with  the  PC-PU-PC  structure. 
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Fig.  16.  Stress  on  wall  for  the  different  cases  such  as  PU,  PC,  PU/PC,  PC/PU,  PC/PU/PC  and 
PU/PC/PU 
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Fig.  17.  Kinetic  energy  on  wall  for  the  different  cases  such  as  PU,  PC,  PU/PC,  PC/PU, 
PC/PU/PC  and  PU/PC/PU 
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